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We apply symmetry considerations to study the possible superconducting order parameters in 
ferromagnetic ZrZn2. We predict that the presence and the location of the superconducting gap 
nodes depend on the direction of magnetization M . In particular, if M is directed along the z 
axis, then the order parameter should always have zeros. We also discuss how to determine the gap 
symmetry in ZrZn2 using ultrasound attenuation measurements. 
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Recently, superconductivity has been found in ZrZn2 
jjj, which a weak itinerant ferromagnet. The most sur- 
prising fact is that the superconductivity occurs only 
in the ferromagnetic phase. The exchange splitting of 
the Fermi surfaces makes a conventional singlet BCS-like 
pairing strongly suppressed. A number of theories have 
been proposed which show how the exchange by spin 
fluctuations can lead to a triplet Cooper pairing both 
in the paramagnetic and the ferromagnetic phases 0, 
or to the enhancement of the superconducting critical 
temperature T c on the ferromagnetic side Another 
feature of the phase diagram is that T c grows as pres- 
sure moves away from the ferromagnetic quantum criti- 
cal point, which can be explained by the exchange- type 
interaction of the magnetic moments of the Cooper pairs 
with the magnetization density [0 . 

Symmetry considerations can identify the possible or- 
der parameters, even in the absence of a firmly estab- 
lished microscopic mechanism of pairing, which is often 
the case for unconventional superconductors. The pres- 
ence of ferromagnetism brings about a number of novel 
features in the symmetry analysis. In this article, we 
give a detailed analysis of the pairing symmetry in ZrZii2 
and discuss its consequences for the superconducting gap 
structure. Our work has some overlap with the recent 
theoretical studies ^] of the gap symmetry in another 
ferromagnetic superconductor UGe2 [Pj. We find that 
the presence and location of the gap zeros depend on the 
direction of magnetization. We also discuss in some de- 
tail the design of ultrasonic attenuation experiments that 
can be used for experimental probing the order parame- 
ter symmetry. 

The symmetry group Q of the system in the normal 
state is defined as a group of transformations which leave 
the system Hamiltonian H$ invariant, i.e. [G,Hq] = 
for all elements G G Q. In non-magnetic supercon- 
ductors, time reversal symmetry K is not broken, and 
Q = S x K x where S is the space group of the 

crystal, and U(l) is the gauge group In contrast, 
in magnetic superconductors, time reversal symmetry is 
broken, and Q = Sm X U{\), where Sm is the magnetic 
space group which is a group of symmetry operations 
leaving both the crystal lattice (the microscopic charge 
density) and the magnetization density M invariant ||. 
For example, if there is a crystal point group rotation R 



which transforms M to —M, then the combined oper- 
ation KR will be an element of Sm, because the time 
reversal restores the original M not affecting the lattice 
symmetry. In the above expressions, it is assumed that 
the space group elements act on the orbital and spin coor- 
dinates simultaneously, which is the case when the spin- 
orbit coupling is present. In the absence of spin-orbit cou- 
pling, the transformations of the orbital and spin spaces 
are independent, so that Q — S or b x G spjrl x 11(1), where 
G S pin = SO(2) is the group of spin rotations about the 
direction of M. 

The crystal structure of ZrZn2 in the absence of ferro- 
magnetism is described by a face-centered cubic Bravais 
lattice, with the Zr atoms forming a diamond structure. 
For cubic ferromagnets, the only two possibilities for the 
easy direction of magnetization are a [001] or a [111] di- 
rection, and these possibilities are analysed in the ar- 
ticle. Since a relatively small magnetic field (0.05T at 
T = 1.75 A') is required to line up the magnetic moments 
along a given direction jD, it is expected that experi- 
ments could be carried out for M parallel to either [001] 
or [111]. The change in the superconducting gap struc- 
ture when M is rotated by an external magnetic field is 
one of the interesting and unusual properties of ZrZn 2 
that could be investigated experimentally. 

If Af is along [001], then Sm is generated by: (i) 
the lattice translations by the primitive vectors of the 
fee lattice: t x = (ffl/2)(l, 1, 0), t 2 = (a/2)(0, 1, 1), and 
*3 = (a/2)(l, 0, 1), where a is the lattice constant; (ii) the 
rotations C^ z about the z axis by an angle 7r/2 followed 
by a fractional translation by a vector r = (a/4)(l, 1,1); 
(iii) the combined rotations KC 2x about the x axis by an 
angle 7r accompanied by the time reversal; and (iv) the in- 
version A The point symmetry of the crystal is described 
by the magnetic group D^C^) = D 4 (C 4 ) x Cj, where 
Cj = {E,I}. The subgroup in parentheses (the unitary 
subgroup) incorporates all symmetry elements which are 
not multiplied by the anti-unitary and anti-linear opera- 
tion KC 2x , i-e. D 4 (C 4 ) = C 4 + KC 2x x C 4 . 

If M is along [111], then Sm is generated by: (i) 
the primitive lattice translations, (ii) the rotations C^ xyz 
about the [111] direction by an angle 2ir/3; (iii) the 
combined rotations KC 2xy about the [110] direction by 
an angle ir accompanied by the time reversal; and (iv) 
the inversion I. The point symmetry is described by 
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the magnetic group D 3 d(C 3i ) = D 3 (C 3 ) x Cj, where 
D 3 (C 3 ) = C 3 + KC 2xy x C 3 . 

Using the standard notation for the space group op- 
erations which combine rotations R and translations t: 
r — > (R\t)r = Rr + t, the transformation rules for the 
spinor wave functions can be written as 

(R\t)Mr) = [D^WU^iR-^r -*)). (1) 

Here a =f, \ is the spin projection on the direction of 
M, and D^ 1 / 2 \R) is the spinor representation of rota- 
tions: for a rotation by an angle 9 around some axis n, 
D(im(R) = uk ] (0) = exp(-i(0/2)(tr-n)). We will also 
need the transformation rules under the time reversal op- 
eration: 



and the inversion: 



Iipa(r) = t/v(-r). 



(2) 



(3) 



In the presence of the exchange field and spin-orbit 
coupling, the single-particle wavefunctions are linear 
combinations of the eigenstates of the spin operator s z : 
(r\tf)) = u(r)\ |) + v(r)\ |). Because the normal state 
Hamiltonian H$ is invariant with respect to the crys- 
tal lattice translations, the eigenfunctions are the Bloch 
waves tjjk(f) corresponding to wave vectors k in the Bril- 
louin zone. If the energy spectrum consists of a single 
band which is doubly degenerate in zero exchange field 
due to the Kramers theorem, then diagonalization of the 
Hamiltonian in the presence of the exchange field results 
in two non-degenerate energy bands e±(k). The corre- 
sponding single-particle wave functions have the form 



(r\k,±) 



,(±: 



(r)| T> 



'(r)| !>• 



(4) 



These states are referred to as the pseudospin states. The 
operations from the group Q conserve the pseudospin in 
the following sense: G\k,±) — exp[i0±(fc, G)]|G or &fc, ±), 
with G or b describing the "orbital" part of the symmetry 
operation, e.g. rotations or reflections, and the unde- 
termined phase factors <j>± coming from the freedom in 
choosing the overall phases of \k, ±) at every point of the 
Brillouin zone. Here we adopt a convention introduced 
in Ref. [|l0| , according to which the pseudospin states (Q) 
transform similar to the spin eigenstates |fe, |) and |fc, j). 

Let us first consider the case M || [001]. From Eqs. 
(Q,pH), the transformation rules for the creation opera- 
tors of electrons in the states (0) are: 



(G 4z |r) 



c k,± 



-i(C iz k)-T =pi7r/4„t 



(KC 2x \0): A4 i± - ±a*cl c ^ fe)± 



(5) 



u fc,± 



-fc,±- 



Here A is an arbitrary c- number. 

The pseudospin states can be used as a basis for con- 
structing the Hamiltonian which takes into account the 



Cooper pairing between electrons with opposite momenta 
k and —k. Treating the Cooper interaction in the mean- 
field approximation, we obtain H — Hq + H SCl with the 
non-interacting part 



Ck+ + e_(fc)ct c fe _ 



(6) 



describing two separate sheets of the Fermi surface cor- 
responding to different pseudospin indices, and 



F - = E E [ A «/3(fc) c L c+ -M +h - c - 

" c*,/3=± 



(7) 



Here A ++ (fc) and A__(fc) represent the superconduct- 
ing order parameters at the "+" and "— " sheets of the 

Fermi surface respectively, and A H (k) = —A |_(— k) is 

the order parameter composed of quasiparticles on differ- 
ent sheets. From the Pauli exclusion principle, A ++ (fc) 

and A__(fc) are odd functions of k, but A H (k) does 

not have a definite parity. Separating the odd and 
the even parts, the order parameter matrix can also 
be cast in a more familiar form A(fe) = (iercr 2 )d(fc) + 
(i<72)ip{k), where d and ip are the pseudospin-triplet and 
the pseudospin-singlet components respectively B. The 
fact that the Fermi surface of ZrZii2 consists of several 
sheets of different topology jllj , does not change our re- 
sults. 

From Eqs. (^|), the band spectra e±(fe) are invariant 
under the operations from the point group D^. Also, we 
obtain the transformation rules for the order parameters 
under rotations C± z : 



A ++ (fc) - -«A ++ (C£ l fc) 
A__(fc) -» +iA__(G 4z 1 fc) 



(8) 



(note the cancellation of the x-dependent phase factors 
on the right-hand side of these equations), and under the 
combined time reversal and rotations KC 2x : 



A ++ (fc) -» A* ++ (C^k) 
.(fc)-A*_(C£fc) 



A_. 

A+- 



(9) 



K-(-c^k). 



In the presence of the exchange band splitting E ex , the 
low-frequency part of the spectrum of excitations (e.g. 
spin fluctuations) responsible for the interband Cooper 
pairing is cut out 0. Since E ex is by far the largest en- 
ergy scale in the system: E ex ~ 5mRy ~ 800K [12], the 
pairing interactions c k ,c_ k _Cfc<_c_fc'. + , which are re- 
sponsible for A^ , are negligibly small jTj| (some of the 

consequences of taking these interactions into account 
will be discussed below). On the other hand, the inter- 
band pairing terms c\. + cL k , Ch'-C-k',- can induce order 
parameters of the same symmetry on both sheets of the 
Fermi surface. We expect the effect of these terms to 
be small at small spin-orbit coupling, because they are 
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TABLE I: The character table and the examples of the odd 
basis functions for the irreducible co-representations of the 
magnetic point group D^Ci). The overall phases of the basis 
functions are chosen so that KCzxfvik) — fr{k). \i,2 are 
arbitrary real constants. 



r 


E 




Hk) 


A 


1 


1 


k z 


B 


1 


-1 


k z [\i(ky + ik x ) 2 + A 2 (fc H — ik x ) 2 ] 


X E 


1 


i 


ky ~\~ ikx 


2 E 


1 


—i 


ky ikx 



absent at zero spin-orbit coupling due to the spin conser- 
vation. 

The superconducting order parameter which emerges 
at T c transforms according to one of the irreducible 
representations T of the normal state symmetry group 
Q. It can be represented as the expansion Ap(fc) = 
J2iVr,ifi{k), where i labels the orbital basis functions, 
and ?7r,i are the order parameter components which en- 
ter, e.g., the Ginzburg-Landau free energy. In our case, Q 
contains the anti-unitary and anti- linear operation KC2 X , 
and, instead of usual representations, one should use 
co-representations of the magnetic point group D4(C4), 
which can be derived from one-dimensional representa- 
tions of the unitary subgroup C4 [Q. The results for 
odd co-representations are listed in Table Q. Note that 
the action of the unitary and anti-unitary orbital sym- 
metry elements on scalar functions f(k) is defined as: 
Rf{k) = /(R^k), and KRf(k) = f^-R^k). 

If the superconductivity appears on the "+" -sheet, 
then the order parameter is A ++ (fc) [in terms of the vec- 
tor order parameter d(k) — d z (k)z + (d+(k)(x — iy) + 
d-(k)(x + iy))/2, it corresponds to d- = d x — id y ]. Using 
Eqs. (|J) and Table EL we obtain the following expressions: 

A ++ . j4 (fc) = iVA^E(k) 
A ++ ,s(fc) = iVBh E {k) 
A ++ i B (k) = ir]i E f B (k) 
A++,2 £; (fc) = ii]2 E f A (k). 



(10) 



The appearance of different representations on the left- 
hand and the right-hand sides of these expressions can be 
easily understood if to look at Eqs. (g). The transformed 
order parameter A ++ has an extra factor —i which comes 
from the rotation of spin coordinates. In terms of d(k), 
this factor is the result of the rotation of the basis spin 
vector x + iy, which transforms according to the 2 E rep- 
resentation. Thus, for instance, the first line of Eqs. jlO| ) 
follows from the fact that A — 2 E x 1 E. So far, we have 
discussed the transformation properties of the order pa- 
rameters (|Io| ) under the rotations from the unitary sub- 
group C4. Because of our choice of the overall phase of 
the basis functions (see the caption to Table pj) and the 
factors i on the right-hand sides of Eqs. (jl(]), the effect 
of the anti-unitary operation KC^x on rjr is equivalent 
to complex conjugation: KC2 X r\v — Vr ■ 



As seen from Eqs. (10) and Table SJ, the order parame- 
ters Aa and A E vanish at the poles of the Fermi surface 
k x = k y = 0, while the order parameters Ai E and A2 E 
vanish at the equator k z = 0. One can prove that these 
gap zeros are not artifacts of our choice of the basis func- 
tions but are imposed by symmetry. Indeed, one of the 
elements of the unitary component C^h of the magnetic 
point group is the basal plane reflection <ih — C2Z x I- 
Therefore, 

ChfA,B(k) = fA,B(k x , ky, -k z ) = - fA.B(k), (11) 

so that fA,B(k x ,k y ,0) = 0, and Ai E (k x ,k y ,0) = 
A2 E (k x ,k y ,0) = 0. Similarly, under a four-fold rotation 
around the z axis: 

C± z fiE, 2 E(k) = fi E 2 E (k y , -k x , k z ) = ±ifi E 2 E (k), 

hence /ib,2b(0, 0, k z ) = 0, and A / i(0,0, k z ) — 
A B (0,0,fc z j = 0. It also follows from Eq. (|ll|) that 
$A{k) and f E {k) go to zero at k z — ±n/a, i.e. at 
the surface of the Brillouin zone, because (k x , k y , ■n/a) 
and (k x ,k v , —n/a) are equivalent points. In order to 
take into account the crystal periodicity leading to the 
presence of these additional gap zeros, one has to rep- 
resent the basis functions as the lattice Fourier series 
/(&) = Sn/™ elfc _R "' wn ere summation goes over the 
sites R n of the Bravais lattice of the crystal. The ex- 
pansion appropriate for an odd order parameter has the 
form 



/(*) = 



sin k ■ R n , 



(12) 



where R n are the sites of a fee cubic lattice, 
which cannot be transformed one into another by 
inversion. In the nearest-neighbor approximation, 
we choose the following set of i?„'s: {Rn} — 
o/2{(101), (101), (Oil), (011), (110), (110)}. Using the 
representation characters from Table SJ, we obtain the ba- 
sis functions which have symmetry-imposed zeros at the 
surface of the Brillouin zone: 



f A (k) 



sin ■ 



k z a 



cos ■ 



f B (k) = sm^y 



cos 



2 



cos 



cos 



kyOb 



fi E {k) = cos — I sin -|- 



+Ai 



iTr / k x a 


k v a\ 




e 4 sin 1 — - — 


— e 




~) 





k z a 



ky (X 



f2 E (k) = cos — sin i sin — 



+A 2 



4 sm 



2 



e 4 sm 



2 



2 



Here X12 are arbitrary real constants. The polynomial 
expressions for the basis functions from Table m are re- 
covered in the limit of a "small" Fermi surface k — > 
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[note that /s(fc) from Table | can be obtained by includ- 
ing the next-nearest-neighbors in the expansion It 
should be noted that these nearest-neighbor results give 
also gap zeros not required by symmetry, e.g. /s(fc) = 
on the plane k x = k y . These "accidental" zeros will be 
removed if higher-neighbor terms are included, but if the 
nearest- neighbor terms turn out to be dominant, experi- 
ment could find indications of these accidental zeros. 

The order parameter on the pseudospin-down sheet of 
the Fermi surface is A__(fc) [in terms of d(k), it cor- 
responds to eL|_ = d x + id y , and the relevant basis spin 
vector is x — iy, which transforms according to the 1 E 
representation] . Its fe-dependence is given by the follow- 
ing expressions: 



A — . A (k) = iriA^E{k) 
A — , B (k) = ir] B f^ E {k) 
A — .i£(fc) = irji E fA(k) 
A — 2 E (k) = iip E f B (k). 



(13) 



If we take into account the interband pairing interac- 
tion of the form c k+ c^_ k + Ck' -C_fe' then both A ++ and 

A are non-zero and correspond to the same irreducible 

co-representation of the magnetic point group. Compar- 
ing Eqs. (10) and Jl5|), we see that, although the orbital 

symmetries of A ++ and A are different, they have 

the symmetry-imposed gap nodes at the same locations 
on both sheets of the Fermi surface. For example, if the 
order parameter corresponds to the A co-representation, 
then both A ++iJ 4(fc) ~ fi E (k) and A__ i x(fc) ~ f 2 E{k) 
have point nodes at k x — k y = 0. 

The gap nodes disappear only if the interband pair- 
fc _c/c_c_fc' !+ are taken into ac- 
These terms induce a non-zero order parame- 



ing interactions c' k+ c 



count. 

ter A_| , whose momentum dependence in the triplet 

channel, according to Eqs. (JsJ) and Table |[ is given by 
A+_ r (fe) ~ fr{k), where F = A,B, X E, or 2 E. To sec 
explicitly how the structure of the nodes in the different 
components of the gap function is translated into zeros 
of the spectrum of Bogoliubov quasiparticle excitations, 
it is necessary to diagonalize the Hamiltonian of Eqs. @ 
and (0). This gives the following condition for the energy 
E(k) of a quasiparticle to be zero at some k: 



elel 



|A_ 



-(|A H 



IA, 



I A H 

|2\ 



|det A| 2 = 0. (14) 



The condition for zeros in the excitation energy on the 
"+" -sheet of the Fermi surface (i.e. at e+(fe) = 0) is that 

A ++ (fc) = A _| (fe) = 0, while for zeros in the excitation 

energy on the "—"-sheet (i.e. at e_(fe) = 0) we must 

have A (fe) = A_| (fc) = 0. Thus a nonzero A_| will 

remove the nodes in the spectrum of elementary excita- 
tions. For example, A H ,A(k) ~ f A (k) does not vanish 

at k x = k y = 0, so that these point nodes should be filled. 
However, as discussed above, this effect is expected to be 
negligibly small. 

In the absence of a complete understanding of the mi- 
croscopic mechanism of the superconductivity in ZrZii2, 



TABLE II: The character table and the examples of the odd 
basis functions for the irreducible co-representations of the 
magnetic point group D3(Cs), lj = e 27rj//3 . The overall phases 
of the basis functions are chosen so that KC2 xy fr(k) = fr(k). 
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C^xyz 




Hk) 
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1 




l E 


1 




(jj* 


— i7r/3r, jL , p ^/3jL 


2 E 


1 


UJ* 


Id 


p^l^h h -1- p~™l Z h 

o I\jx — ivy \^ C /1j2 



one cannot tell which order parameter from the lists (|Io| ) 
and ( p^ ) corresponds to the highest critical tempera- 
ture. For example, if the superconductivity is due to the 
exchange by spin fluctuations, then, at vanishing spin- 
orbit coupling, the order parameters A ++ ^a, A ++ , Bl and 
A_)__|_ 2 E correspond to the p-wave equal-spin-pairing su- 
perconducting states studied in Ref. in these terms, 
A ++ ^i E corresponds to /-wave pairing. In Ref. a 
simple phenomenological model of the phase diagram of 
ZrZn2 was proposed. The basic idea was that the under- 
lying order parameter is a vector quantity transforming 
according to a three-dimensional representation of the 
cubic group. Then, the exchange-type interaction of the 
magnetic moments of Cooper pairs with the ferromag- 
netic magnetization splits the superconducting critical 
temperature and lowers the dimensionality of the order 
parameter from three to one. In this model, the order pa- 
rameters Ai E and A2 E are the possible ones, and the ex- 
perimental determination of the gap symmetry will thus 
be helpful in assessing the validity of the model. 

A similar analysis can be done if the easy axis for mag- 
netization is [111]. In this case, the band spectra are 
invariant under the operations from the group D 3t ;, the 
relevant magnetic point group is D 3 (C 3 ), and the trans- 
formation rules (j^) and (Q) for the order parameter are 
replaced by 



A ++ (fc)^ e - 2 -/ 3 A ++ (C7 3 - 1 yz fc) 
A__(fc)^ e +2 -/ 3 A__(C7 3 -4fc) 
A+_(fc)^A+_((7 3 -^fe), 



(15) 



under rotations C 3a; j /2 , and 



A ++ (fe) - A* ++ (C^ y k) 
A__(fc)-A*_(C£„fc) 
A + _(fc)^A;_(-C 2 ^fc), 



(16) 



under the combined operation KC^xy 

Using Table |fi|, we obtain the following fc-dependences 
of the order parameter at the pseudospin-up sheet: 



A + +, A (k) = iVAfi E (k) 
A ++ ,iE(fe) = irii E f2 E (k) 
A ++ 2 E (k) = ir/2 E f A (k). 



(17) 
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As above, the factors i here guarantee that the action of 
KC2xy on ?7r is equivalent to complex conjugation. 

A consequence of the above results is that there are 
no gap nodes for A2 E , whereas A^ and Ai E have point 
nodes where the line k x = k y = k z cuts the Fermi surface. 
For completeness, we also give the expressions for the 
basis functions of the magnetic point group 03(03) in 
terms of the lattice Fourier series in the nearest-neighbor 
approximation: 

f A (k) = S+ + S+ + S+ 

(Sf + Sa + S3) 
fi E (k) = lu*S+ +ojS} + 

+1X2 (to* Si + + S3) 
f2 E (k) = ujS? +uj*S} + £$ 

+1X3 + u>*S2 + S3) , 

where Sf' = sin(fc x a/2 ± k y a/2), — shi(k y a/2 ± 
k z a/2), S% = sin(fc z a/2 ± k x a/2), and Ai,2,3 are arbi- 
trary real constants. 

For the order parameter at the "—"-sheet of the Fermi 
surface, 



A — lA (k) = irj A f2 E (k) 
A — x E (k) = i^EfAik) 
A — 2 E (k) = irj2 E fi E (k). 



(18) 



There are no gap nodes for Ai E , but A a and A2 E have 
point nodes at k x — k y = k z . If the interband hybridiza- 
tion of the form c' k , c_ fe , Ck'-C-k',- is taken into ac- 
count, the order parameters are non-zero on both sheets 
of the Fermi surface. From Eqs. (|l7]) and (|l8|), we see 
that both A ++ , A (fc) - /i B (fe) and A__, A (fc) - f2 E (k) 
vanish on the line k x = k y = k z . The gap nodes dis- 
appear only in the presence of the interband pairing 
c\. + c^_ k Cfe'—C—fe/^, which induces the order parameter 

A+_ r (fe) ~ /r(fe), where T = A, X E : or 2 E. Again, we 

expect A^ to be negligibly small in the presence of the 

large exchange field. 

The presence of the gap nodes would manifest them- 
selves in power-law temperature dependences of the ther- 
modynamic and transport properties ||. For example, 
the electronic specific heat at low temperatures should 
be C(T)/T = 70 + 71T for the line nodes, and C(T)/T = 
7o + 7iT 2 for the point nodes. The temperature- 
independent contributions on the right hand side of these 
equations come from the normal excitations at the un- 
paired sheet of the Fermi surface. If the magnitudes of 

both order parameters A ++ are A are comparable (we 

expect this to be the case only if the spin-orbit coupling 
is strong enough) , then 70 is absent and a power-law de- 
pendence should be observed. 

One of the most powerful methods of determining the 
presence and the location on the Fermi surface of gap 
nodes in unconventional superconductors has been ul- 
trasonic attenuation (l|, p, ||, The method 
is based on finding which sound waves are particulary 



weakly attenuated by the nodal quasiparticles. Nodal 
quasiparticles are "inactive" in attenuating a particu- 
lar sound wave if the electron-phonon interaction for 
the nodal quasiparticle with the particular sound wave 
is zero. Symmetry arguments determining the inactive 
nodes have been developed in a previous article jl!| . We 
refer to that article for a detailed treatment of the basic 
ideas, and give here the extension of the arguments which 
is necessary for treating magnetic groups. 

If a symmetry operation of the crystal (i.e. an ele- 
ment of its magnetic point group) leaves the wave vector 
k characterizing a given electron state invariant, then 
the interaction of this electron with certain phonons can 
be shown to be zero. Consider a phonon of wavevector 
q and polarization direction e. The interaction of the 
given electron with the given phonon can be shown to be 
zero if the symmetry operation causes an odd number of 
changes of sign of the two quantities iq and e (the factor i 
is important because the time reversal operation contains 
complex conjugation). The transformation rules for fc, q, 
and e are: (i) under the point-group operations R, they 
transform like polar vectors, i.e. k — > i? _1 fc, etc.; and 
(ii) under the combined operations KR, k — » —R~ 1 k, 
e — > R~ 1 e, and q — > —R~ 1 q, so that iq — > iR _1 q. For 
example, suppose that the magnetic group contains the 
symmetry element KC2 X - The wave vectors fc of elec- 
trons lying in the k x — plane are invariant under KCix ■ 
According to the rule just stated, these electrons have 
zero interaction with transverse phonons having their 
wave vectors along the x axis because, under the opera- 
tion KC2 X , iq remains invariant, but e changes sign. 

As shown above, one should expect the order pa- 
rameter in ZrZn2 to have nodes if M \\ [001], when 
the magnetic point group is 04(04). The gap nodes 
are always active for longitudinal sound waves. If the 
order parameter has point nodes (A^ or A^ in Eqs. 
(pp|)), then these point nodes are inactive for the trans- 
verse waves T100 and T110 polarized either in the basal 
plane or along [001], and also for the waves T001 polar- 
ized either along [100] or [110]. (By definition a Thkl 
sound wave is a transverse wave having its wave vector 
q along the [hkl] direction.) If the order parameter has 
line nodes in the plane k z — (Ai E or A2 E in Eqs. 
(pp|)), then these line nodes are inactive for the trans- 
verse waves T100 and T110 polarized along [001], and 
also for T001 waves polarized either along [100] or [110]. 
Note that the attenuation of the T100 and T110 waves 
polarized in the basal plane can be used to distinguish 
between the equatorial line nodes and the point nodes, 
because the former are active, but the latter are inactive. 
The presence of unpaired electrons on one of the sheets 
of the Fermi surface will not cause difficulties in sym- 
metry determination by ultrasonic attenuation as this 
will simply make a contribution to the low-temperature 
temperature-independent background, which is easily 
distinguished from the temperature-dependent contribu- 
tion of the gapped sheet (or sheets) of the Fermi surface. 

To summarize, we have studied the symmetry of the 
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superconducting order parameter in ZrZii2. If the spin- 
orbit coupling is weak then superconductivity should ap- 
pear on only one of the sheets of the Fermi surface. The 
interband scattering can, in principle, induce non-zero 
order parameters on other sheets and also fill the gap 
zeros, but we expect these effects to be small. The sym- 
metry of the order parameter depends on the direction 
of the easy axis for magnetization. If M || [001], then 
the magnetic point group is D4(C4), and the order pa- 
rameter goes to zero on the line k x — k y — for the gap 
symmetries A and B, or on the planes k z = 0, ±7r/a for 
the symmetries 1 E and 2 E, on both sheets of the Fermi 
surface. The positions of the gap zeros can be probed 
by ultrasonic attenuation measurements, and to assist in 
the design of appropriate experiments we have given a 
detailed discussion of the zeros of the electron-phonon 



interaction in ferromagnetic ZrZn2 which are imposed by 
the magnetic point symmetry. If M || [111], then the 
magnetic point group is D 3 (C 3 ), and the order param- 
eter has point zeros on the line k x — k y — k z on both 
sheets of the Fermi surface, for the gap symmetry A, and 
on one of the sheets, for the symmetries 1 E and 2 E. It 
should be possible to fix the magnetization density M 
along an arbitrary crystallographic direction by the ap- 
plication of an external magnetic field, and hence to de- 
termine the gap structure for ZrZn2 for M along both 
[001] and [111] and to find the changes that occur when 
M is rotated from [001] to [111]. 

We acknowledge a stimulating discussion with Louis 
Taillefer and support from the Canadian Institute for 
Advanced Research and from the Natural Sciences and 
Engineering Research Council of Canada. 
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